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ABSTRACT 

Angular asymmetries have been measured in polarized decays to bbg collected by the 
SLD experiment at the SLC. A high purity bbg event sample is selected by utilizing B 
lifetime information given by the SLD CCD pixel vertex detector and the stable micron-size 
SLC beams, and the b- and b-jets are identified using lifetime information and momentum- 
weighted track charge. The forward-backward asymmetry is observed in the b-quark polar 
angle distribution, and the parity-violation parameter is measured to test the Standard Model. 
Two angular correlations between the three-jet plane and the polarization are studied. 
The CP-even and T-odd, and the CP-odd and T-odd, angular asymmetries are sensitive to 
physics beyond the Standard Model. The latter requires tagging both the b- and b-jet. We 
measure the expectation values of these quantities to be consistent with zero and set limits 
on the correlations at the 5% level. 
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1. Introduction 



The forward-backward polar-angle asymmetry in hadronic decays to two jets has 
been investigated extensively at SLC and LEP to test the predictions of the electroweak 
theory of parity-violation in the Z^qq coupling. In particular, at SLC where the elec- 
tron beam is highly polarized, the left-right-forward-backward asymmetry removes the 
dependence on the Z°e"'"e~ coupling and is directly sensitive to the Z^qq coupling. The 
experimental results are found to be consistent with the theory to within experimental 
uncertainties of a few percent [0]. Hadronic Z^ decays to three jets can be interpreted 
in terms of the fundamental process Z^ —>■ qqg where one of the quarks has radiated a 
gluon. Given the success of the electroweak theory in predicting the two-jet polar-angle 
asymmetry, similar angular asymmetries can be measured in three-jet events to test 
Quantum Chromodynamics (QCD). The Z^ ^bbg final state is particularly interesting 
as a search-ground for possible new physics processes beyond the Standard Model, and 
a high purity sample can be obtained with high efficiency due to the large mass and 
long lifetime of B-hadrons. Here we report the first experimental study of angular 
asymmetries in polarized Z^ decays to bbg. 



2. Angular asymmetries in qqg 

The differential cross section for e"'"e~ qqg can be expressed as 
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where x and x are the scaled momenta of the quark and anti-quark, respectively, 6 is 
the polar angle of the thrust axis 0] w.r.t. the electron beam, and x is the azimuthal 
angle of the event plane w.r.t. the quark-electron plane. Here the thrust axis is 



defined so that it is parallel to the quark direction if the quark has the highest energy, 
and anti-parallel to the anti-quark direction if the anti-quark has the highest energy. 
The cross-section consists of six terms, each of which may be factorized into three 
contributions: 1) event orientation factor in terms of 9 and x; 2) (PcTi/dxdx (i=U,..,A) 
determined by QCD; and 3) hy determined by the fermion electroweak couplings and 
beam polarization. While the first four terms are P-even, the last two terms are P-odd, 
and are sensitive to any parity-violating interactions at the Z^qq and gqq vertices. In 
addition to these six terms, the most general differential cross section can have three 
more terms that are odd under time-reversal 0]. Being T-odd, however, these terms 
vanish at tree level in a theory that respects CPT invariance. 

Recently Burrows and Osland have proposed new QCD tests in terms of the event 
orientation angles |^. Integrating over scaled momenta and the polar angle distri- 
bution of the thrust axis can be expressed as 
d(T 

a(cos 9) = oc (1 - Pe- ■ Ae) (1 + a cos^ 9) + 2Ap{P^- - A^) cos 9, (2) 
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where is the parity violation parameter: 

Ap = ^^Aj. (3) 

Here a = f^j^ff^"? = I '^dE^^^^^ ^"^^ -^e- is the signed electron beam polarization, 
and Ae {Aj) is the electroweak coupling of the to the initial (final) state, given 
by Ai = 2viai/{yf + of) in terms of the vector Vi and axial-vector couplings. By 
manipulating the polarization sign of the electron beam the left-right-forward-backward 
asymmetry, App-, is directly sensitive to the asymmetry parameter Ap, 
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Similarly, by integrating over cos 6*, x and x, the azimuthal- angle distribution can be 
expressed as 

d(j Stt 
27r^ OC (1 - Pe- ■ Ae)(l + /5 COS 2^) - -^A'p^P,- - A,) cosx, (5) 



where A'p is the parity violation parameter: 

A', = -^Aj, (6) 

with (3 = ■ Given the value of the electroweak parameter Af, measurement of the 

angular asymmetry parameters Ap and A'p in ^qqg events allows one to test the 
QCD prediction for apjidu + ol) and (JaI{pu + ^l)- Furthermore, the ratio ApjJ^p 
is independent of Aj and is proportional to ctaIo-p- 

The differential cross-section can also be expressed in terms of the polar angle uj 
of the vector n normal to the event plane w.r.t. the electron beam direction, where 
cosct;=sin6'sinx: 

dcT X6 

OC (1-Pe- ■ Ae){l + COi? U) + —Ar^Pe- -Ae)cOSUJ, (7) 
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with 7 = ~ I aJ+273^L+2/3CTT " ^^^^ ordcr in perturbative QCD, ctl = 2(Tt, yielding 
7 = — |, and A^ = 0. The second term is one of the three T-odd terms mentioned above, 
and appears as a forward-backward asymmetry of the event-plane normal relative to 
the polarization axis. The left-right-forward-backward asymmetry in cos u can also 
be defined by a similar double asymmetry as Eq. 4, and is directly proportional to 
the T-odd parameter At- The vector normal to the event plane can be defined in two 
ways: 1) the three jets are ordered according to their energies, and the two highest 
energy jet momenta are used to define n = pi x and 2) the quark and anti-quark 
momenta are used to define n = Pg x pg. The asymmetry term is CP-even in the 
first definition, and CP-odd in the second. The first definition does not require jet 
fiavor identification, and we have studied the asymmetry for inclusive hadronic Z° 
decays ||^. The second definition requires tagging both quark- and antiquark-jets. In 
both cases, in the Standard Model the asymmetry vanishes identically at tree level, 
but higher-order processes yield non-zero contributions for e"'"e~ -^bbg. However, due 
to various cancellations, these contributions are found to be very small at the Z^ 
resonance and yield values of the asymmetry parameter \At\ < 10~^ 0. Measurement 
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of the asymmetry in cosu; is hence potentially sensitive to physics processes beyond 
the Standard Model 



3. Event and Track Selection 

The SLAC Linear Collider (SLC) collides longitudinally polarized electrons with unpo- 
larized positrons at a center-of-mass energy of 91.2 GeV. The electron polarization di- 
rection is randomly reversed pulse-by-pulse, reducing systematic effects on polarization- 
dependent asymmetries. The magnitude of the average electron-beam polarization was 
0.63 in 1993, 0.77 in 1994-1996, and 0.73 in 1997-1998. 

The measurement was performed with the SLC Large Detector (SLD) using approx- 
imately 550,000 decays collected between 1993 and 1998. A general description of 
the SLD can be found elsewhere [^. Charged particle tracking and momentum anal- 
ysis is provided by the central drift chamber (CDC) |10| and the CCD-based vertex 
detector (VXD) |ll| in a uniform axial magnetic field of 0.6 T. About 70% of the 
data were taken with a new vertex detector (VXD3) installed in 1996, and the rest 
with the previous detector, VXD2. Particle energies are measured in the liquid argon 
calorimeter (LAC) [|1^ and in the warm iron calorimeter . 



In the present analysis the hadronic event selection, three-jet reconstruction, and 
b-tagging were based on charged tracks. A set of cuts was applied to the data to select 
well-measured tracks and events well contained within the detector acceptance fLi 



Events were required to have (i) at least 7 charged tracks; (ii) a visible charged energy 
of at least 20 GeV; and (iii) a thrust axis polar angle satisfying |cos6'j'| <0.71, 
which was reconstructed using the LAC. Charged tracks reconstructed in the CDC 
were linked with pixel clusters in the VXD by extrapolating each track and selecting 
the best set of associated clusters. The average efficiency of reconstruction in the 
CDC and linking to the correct set of VXD hits is 95% (94%) for the region | cos^^| < 
0.85 (0.74) [ITB]. The momentum resolution of the combined CDC and VXD systems 
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is = (.01)^ + (.0026p_l)^, where p± is the transverse momentum in GeV/c 

w.r.t. the beamhne. 

The centroid of the micron-size SLC Interaction Point (IP) in the r0 plane is re- 
constructed with a measured precision of a^.^^ ~ 5/im (7/im) using tracks in sets of ~30 
sequential hadronic decays. The z position of the IP is determined on an event- 
by-event basis with a precision of crf^ ~32/im (52/im) using the median z position of 
tracks at their point-of-closest approach to the IP in the r0 plane. The track impact 
parameter resolution at high momentum is ll/im (ll/xm) in the plane perpendicular 
to the beam axis (r0 plane) and 23/im (38/im) in the plane containing the beam axis 
{rz plane). 

A set of "quality" tracks for use in heavy quark tagging was selected. Tracks 
measured in the CDC were required to have >40 hits, with the first hit at a radius r < 
39 cm, a transverse momentum p±_ > 0.4 GeV/c, a good fit quality {x^/Ndof < 5), 
and to extrapolate to the IP within 1 cm in r0 and 1.5 cm in z. Tracks were required to 
have at least one associated VXD hit, and a combined CDC- VXD fit with x^/^dof < 
5. Tracks with an r0 impact parameter 6 > 3 mm or with an impact parameter error 
as > 250 fim were removed. Tracks from identified 7 conversions and or A° decays 
were also removed. 



4. bbg Analysis 

Three-jet events were selected and the three momentum vectors of the jets were recon- 
structed. Although the parton momenta are not directly measurable, at ^/s ^91 GeV 
the partons usually appear as well-collimated jets of hadrons. Jets were reconstructed 
using the "Durham" jet algorithm [|r^. Planar three-jet events were selected by requir- 
ing exactly three reconstructed jets to be found with a jet-resolution parameter value 
of ?/c=0.005, the sum of the angles between the three jets to be greater than 358°, and 
that each jet contain at least two charged tracks. A total of 75,000 events satisfied 
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these criteria. 

Such jet algorithms accurately reconstruct the parton directions but measure the 
parton energies poorly fl^. Therefore, the jet energies were calculated by using the 



measured jet directions and solving the three-body kinematics assuming massless jets, 
and were then used to label the jets such that Ei > E2 > E3. The energy of jet 1, for 
example, is given by 

^ r- sin6'23 , , 

^1 = ^ ■ n M ■ n \ ■ n ^ ^ 
sm 612 + sm 6^23 + sm 6^31 

where 6ki is the angle between jets k and /. 

To select bbg events the long lifetime and large invariant mass of B-hadrons was 
exploited. A topological algorithm was applied to the set of quality tracks in each 
jet to search for a secondary decay vertex. Vertices were required to be separated from 
the IP by at least 1 mm and to contain at least two tracks. Monte Carlo studies show 
that the probability for reconstructing at least one such vertex was ~ 91% (77%) in 
bbg events, ~ 45% (26%) in ccg events, and ~ 2% (2%) in light quark events. Once 
a vertex was found, additional tracks consistent with coming from the vertex were 
attached in an attempt to reconstruct the invariant mass of a B-hadron. A vertex axis 
was formed by a straight line joining the IP and the vertex, which was located at a 
distance D from the IP. For each quality track the distance of closest approach, T, 
and the distance from the IP along the vertex axis to the point of closest approach, 
L, were calculated. Tracks with T < 1 mm, and L/D > 0.25 were attached to the 
secondary vertex, and the vertex invariant mass, Mch, was calculated assuming each 
track was a charged pion. Due to neutral decay products the total momentum vector 
of the tracks and the vertex axis were typically acollinear. To account for the missing 
neutral particles, an additional component of transverse momentum Pt, defined by 
the projection of the total momentum vector perpendicular to the vertex axis, was 



added to yield M = ^M^,^ + P^ + \Pt\ 0. Figures 1(a), (b), and (c) show the 
distributions of this Pr corrected vertex mass for jet 1, 2, and 3, respectively. An event 
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was selected as bbg if at least one jet contained a vertex with M > 1.5 GeV/c^. A 
total of 14,658 events satisfied this requirement and were subjected to further analysis. 
Monte Carlo studies show that this selection is 84% (69%) efficient for identifying a 
sample of bbg events with 84% (87%) purity, and containing 14% (11%) ccg and 2% 
(2%) light-flavor backgrounds. 

The identiflcation of each jet was based on the momentum-weighted jet charge and 
r0 impact parameter techniques. The momentum-weighted charge was calculated for 
each jet: 

Qj = T.^^\Pi-i\''^ (9) 
where k=0.5, i is the unit vector along the event thrust axis, and and Pi are the 
charge and momentum of the i*'* track associated with jet j. We then examined the 
difference in the momentum- weighted jet charge, Qdiff = Qi — Q2 — Qz- If this quantity 
was negative (positive), jet 1 was tagged as the b-jet (b-jet). The jet flavor was tagged 
by counting the number of "signiflcant" tracks with normalized impact parameter w.r.t 
the IP djad > 3. Figures 2(a), (b), and (c) show the distributions of the number of 
signiflcant tracks found in jets 1, 2, and 3, respectively, in the b-tagged events. Jet 1 
was chosen as the gluon-jet only if jet 1 had no signiflcant track and both jet 2 and 
3 had at least one signiflcant track. Jet 2 was chosen as the gluon-jet if jet 2 had 
no signiflcant track and jet 3 had at least one signiflcant track. Otherwise, jet 3 was 
chosen as the gluon-jet. 



5. Monte Carlo Simulation 

A Monte Carlo simulation of hadronic Z° decays combined with a simulation of 
the detector response was used to study the quality of the jet reconstruction, the b- 
tagging efficiency and purity, and the efficiency of the jet flavor identiflcation. The 
JETSET 7.4 PI event generator was used, with parameter values tuned to hadronic 



e"'"e annihilation data pT|, combined with a simulation of B hadron decays tuned p2 
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to T{AS) data and a detector simulation based on GEANT 3.21 |]23[|. For those events 
satisfying the three-jet criteria, exactly three jets were reconstructed at the parton 
level by applying the jet algorithm to the parton momenta. The three parton-level jets 
were associated with the three detector-level jets by choosing the combination that 
minimized the sum of the angular differences between the corresponding jets, and the 
energies and charges of the matching jets were compared. 

For the T-odd asymmetry analyses the vector normal to the jet plane is measured 
in two ways: 1) using the two highest energy jets, and 2) using identified b- and b-jets. 
In the first method, where the jets are labeled according to their energy, six detector- 
jet energy orderings are possible for a given parton-jet energy ordering. For the three 
cases where the energy ordering of any two jets does not agree between parton and 
detector levels, the direction of the jet-plane normal vector is opposite between the 
parton level and detector level and coscj will be measured with the wrong sign. The 
average probability of measuring cos u with the correct sign in this analysis is estimated 
from the simulation to be 76% (76%). In the second method, where both b- and b- 
jets are identified, the gluon-jet must be tagged correctly, and furthermore, the charge 
assignment of the b- and b-jets must be correct. The average probability of identifying 
the gluon-jet correctly is 91% (88%), and combined with the correct-charge assignment 
probability determined by the self-calibration technique described in the next section, 
the average probability of measuring cosu with the right sign is 64% (63%). 



6. Angular Asymmetries 

Figures 3(a) and (b) show the observed cos6' distributions of the signed-thrust axis for 
event samples collected with left- and right-handed electron beam, respectively. The 
histograms show the backgrounds estimated using the simulation. The thrust axis is 
signed so that it points towards the hemisphere containing the b-tagged jet. The cos 6' 
distribution may be described by 
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^008 6* 

2 (P,- - A,) cose [Ap,h h (2/°'-^'=^*'^ - 1) + 

Ap,, /, (2p— - 1) + Ap,„,, (1 - A - /e) (2p'=°^'^'^^*'"'^^ - 1) ], (10) 

where /c, fuds are the fractions of hbg, ccg, and hght quark events in the sample, 
determined from the Monte Carlo simulation, and p'^°^"^<^'^*^^ ^ pCorrect,c^ ^correctMs ^j.^ 
the probabilities to tag the parton charge correctly for hbg, ccg, and the light quark 
events, respectively. The correct-charge probability is calculated as a function of the 
measured jet charge difference \QdiSf \ as p™^''^^* = + e""'*^*^^'). The quantity a 
is a parametrization of how well the momentum-weighted charge technique signs the 
thrust axis direction. While o;, and auds for ccg and light-quark backgrounds were 
calculated from the simulation, for bbg events was determined from data using a 
self-calibration technique |2^. Using the measured widths, cr^i// and agumi of the 



Qdiff (= Q1-Q2- Q3) and Qsum (= Q1 + Q2 + Q3) distributions: 



2>lff-(l + A)V2 



'diff K-^ ~ '-'sum , , 

(1 + V^(^sum 

where the hemisphere correlation A = 0.027 was calculated from the simulation. This 
yielded at = 0.218±0.021 (0.255±0.032) averaged over cos 6*. On average the correct- 
charge assignment probability for bbg events is 68% (67%). The asymmetry parameters 
Ap^c and Ap^u^s (Eq. |^) for charm and light-quark backgrounds were calculated from 
the simulation based on the Standard Model. A maximum-likelihood fit of Eq. |T0| is 
performed to extract Ap^b- We found 

Ap^b = 0.847 ± 0.049, (PRELIMINARY) (12) 

where the error is statistical only. Assuming the Standard Model expectation of Ab = 
0.94 for sin^ 6*10=0.23, the measured value of Ap^b yields 

, ^ = 0.906 ± 0.052(stat.). 
10 



This value is consistent with the O^a"^) QCD expectation of crp/{cru + ctl) = 0.93, 
calculated using the JETSET 7.4 event generator . 



Figures 4(a) and (b) show the x distributions for event samples collected with left- 
and right-handed electron beam, respectively. The x distribution may be described by 



^ = (l-P,--A,)(l+/5cos2x) - 
dX 



{Pe--A,) COSX [A'p^JbP'AP + 



A'pJcPIp + A'p,u,s (1 - A - /c) Pap' ], (13) 

where Pap, Pap, P^j^ are the analyzing powers for bbg, ccg, and light quark events, 
respectively, and are function of the probability to tag the parton charge correctly, Pchg, 
and the probability to tag the gluon-jet correctly, Pgiu, given by 



Pap — PchgPglu ^ (1 ~ Pchg)Pglu — Pchg{^ — Pglu) + (1 ~ Pchg){.^ — Pglu) ■ (14) 

A maximum-likelihood fit of Eq. [13| is performed to extract v4p^. We found 

A'p^^ = -0.013 ± 0.033, (PRELIMINARY) (15) 

where the error is statistical only. Assuming the Standard Model expectation of Aj, = 
0.94 for sin^ 6^^=0.23, the measured value of yields 

— — — = -0.014 ± 0.035(stal). 

(Tu + ctl 

This value is consistent with the 0{al) QCD expectation of ca/ {o'u + o^l) = —0.064, 
calculated using the JETSET 7.4 event generator . 



Figures 5(a) and (b) show the left-right-forward-backward asymmetry of the cosu 
distribution for the two definitions: (a) pi x p^, and (b) pf, x pj. No asymmetry is 
apparent. The costu distribution may be described, assuming no asjTumetries in the 
ccg and light-quark backgrounds, by 
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— ^ oc (1 - Pe- ■ v4e)(l - -cos'c^) + — [P^- " A^) At fb Pap cosu, (16) 
a cos a; 3 9 

where fb is the fraction of bbg events in the sample, and the analyzing power, Pap, 
represents the probability of correctly signing the vector normal to the event plane. 
In the first case Pap is given by the probability of correct energy-ordering. Pap = 
j^2pCorrect_]^-j^ and in the second case it is the probability of correct-sign assignment com- 
bined with the tagging efficiency of the gluon-jet. We performed maximum-likelihood 
fits of Eq. |l6]to the cosu distributions to extract the parameters A^, for the CP-even 
case, and A^, for the CP-odd case. 
We found 

A^ = -0.012 ± 0.013, (PRELIMINARY) 

= -0.033 ± 0.023, (PRELIMINARY) 

where the error is statistical only. In both cases the T-odd contribution is consistent 
with zero within the statistical error and we calculate limits of 

-0.038 <A^< 0.014 @ 95% C.L., (PRELIMINARY) 

-0.077 <A^< 0.011 @ 95% C.L.. (PRELIMINARY) 
The results of these fits are shown in Figures 5(a) and 5(b). 

7. Systematic Errors 

Table |l| summarizes the systematic errors on the forward-backward asymmetry analysis 
of the signed thrust-axis. The largest systematic error was due to the statistical uncer- 
tainty in the ab determination using the self- calibration technique. This error would 
decrease with a larger data sample. The systematic error in the inter-hemisphere 
correlation A was due to the limited statistics of the Monte Carlo simulation. The 
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systematic error in the tag composition was due to the heavy quark physics model- 
ing. In bb events we have considered the uncertainties on: the branching fraction for 
-^bb, the B hadron fragmentation function, the rates of production of B"^, B^ and 
Bg mesons, and B baryons, the hfetimes of B mesons and baryons, and the average 
B hadron decay charge multiphcity. In cc events we have considered the uncertainties 
on: the branching fraction for Z'^ -^cc, the charmed hadron fragmentation function, 
the rates of production of D^, D'^ and Dg mesons, and charmed baryons, and the 
charged multiphcity of charmed hadron decays. We have also considered the rate of 
production of ss in the jet fragmentation process, and the production of secondary 
bb and cc from gluon sphtting. The systematic error in the detector modeling results 
from discrepancies between data and Monte Carlo in tracking efficiency and resolution. 
The systematic errors on A'p and At are negligibly small as the uncertainty diminishes 
with the asymmetry itself. 

8. Conclusions 

In conclusion, we have made the first angular asymmetry measurements in polarized 
Z^ decays to bbg. From the forward-backward polar angle asymmetry of the signed- 
thrust axis we have measured the parity violation parameter Ap = 0.847 ± 0.049 (stat.) 
± 0.060 (syst.). Prom the azimuthal angle asymmetry we have measured the second 
parity violation parameter A'p — —0.013 ± 0.033 (stat.) ± 0.002 (syst.). Assuming the 
Standard Model expectation of Af, — 0.94, the QCD factors for bbg events are measured 
to be ap/{au + aL) = 0.906 ±0.052 (stat.) ± 0.064 (syst.), and (Ta/(<^c/ + (Xl) = -0.014 
±0.035 (stat.) ± 0.002 (syst.), which are consistent with the 0{a'^) QCD expectations. 
We find the T-odd asymmetry to be consistent with zero, and we set 95% C.L. limits 
on the asymmetry parameter —0.038 < ^4^ < 0.014 for the CP-even case and —0.077 
< A^ < 0.011 for the CP-odd case. All results are preliminary. 
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Table 1: Contributions to the relative systematic error on Ap. 



Error Source 


SAp/Ap 


ab 


5.7% 


Monte Carlo statistics on A 


1.0% 


Tag Composition 


3.9% 


Detector Modeling 


1.0% 


Beam Polarization 


0.8% 


Total 


7.1% 



Figure Captions 

Figure 1. PfCorrected vertex mass distribution for (a) highest energy jets, (b) second- 
highest energy jets, and (c) lowest energy jets. The histograms are Monte Carlo simu- 
lations; the flavor compositions of the simulations are indicated. 

Figure 2. Numbers of significant tracks in b-tagged events for (a) highest energy jets, 
(b) second-highest energy jets, and (c) lowest energy jets. The histograms are Monte 
Carlo simulations; the flavor compositions of the simulations are indicated. 

Figure 3. Polar-angle distribution of the signed-thrust axis direction with respect to 
the electron-beam direction for (a) left-handed and (b) right-handed electron beam. 
The histograms are Monte Carlo estimations of the backgrounds. 

Figure 4. Azimuthal- angle distribution of the signed-thrust axis direction with respect 
to the electron-beam direction for (a) left-handed and (b) right-handed electron beam. 
The histograms are Monte Carlo estimations of the backgrounds. 
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Figure 5. Left-right-forward-backward asymmetry in polar-angle distribution of the 
vector normal to the event plane for (a) CP-even case, and (b) CP-odd case. The sohd 
curve is the best fit to the data sample, and the dashed curves correspond to the 95% 
C.L. limits. 
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